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, Direct phonon spin-lattice relaxation of an electron qubit bound by a donor impurity or quantum 

' dot in SiGe heterostructures is investigated. The aim is to evaluate the importance of decoherence 

' from this mechanism in several important solid-state quantum computer designs operating at low 

temperatures. We calculate the relaxation rate 1/Ti as a function of [100] uniaxial strain, tem- 
(~| ■ perature, magnetic field, and silicon/germanium content for Si:P bound electrons. The quantum 

^ ' dot potential is much smoother, leading to smaller splittings of the valley degeneracies. We have 

estimated these splittings in order to obtain upper bounds for the relaxation rate. In general, we 
. find that the relaxation rate is strongly decreased by uniaxial compressive strain in a SiGe-Si-SiGe 

' quantum well, making this strain an important positive design feature. Ge in high concentra- 

tions (particularly over 85%) increases the rate, making Si-rich materials preferable. We conclude 
that SiGe bound electron qubits must meet certain conditions to minimize decoherence but that 
^ ' spin-phonon relaxation does not rule out the solid-state implementation of error-tolerant quantum 

C/3 ' computing. 

PACS numbers: 3.67.Lx,85.35.Gv,72.25.Rb 
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I. INTRODUCTION 



The prospect of quantum computing (QC) has caused great excitement in condensed matter physics. If a set of 
'"^ , qubits can be maintained in a coherent, controllable many-body state, certain very difficult computaiional problems 
^ ■ become tjcactable. In particular, successful QC would mean a revolution in the areas of cryptographyEJ and data-base 
O I searching3. In addition, it would mean a great advance in general technical capabilities, since the control of individual 
I ^ I ■ quantum systems and their interactions would represent a new era in nanotechnology. 

However, from a practical point of view, a dilemma presents itself immediately. On the one hand, one wishes to 
control quantum degrees of freedom using external influences, since that is how a quantum algorithm is implemented, 
^ ' and to measure them, since that is the output step. On the other hand, the system must be isolated from the 
environment, since random perturbations will destroy the quantum coherence that is the whole advantage of QC. 
This is the isolation-control dilemma, and it leads to a very rough figure-of- merit F for any quantum computer. If 
] we define the decoherence time r as the time it takes to lose quantum coherence, and the clock speed s (roughly the 
' inverse of the time to run a logic gate), then the figure of merit is F ^ st and a practical machine should satisfy 
P\| F > 10** — 10^ at least. If the clock speed is limited only by standard electronics, then we may be able to achieve 
I s w 10® Hz. This would imply that r = 1 ms is a lower limit for the decoherence time. 

The dilemma has not yet been solved, though a number of solutions have been proposed. A particularly attractive 
solution is to use spin degrees of freedom as qubits. Nuclear spins interact relatively weakly with their environment 
because the coupling, proportional to the magnetic moment, is small. Yet there has grown up a sophisticated 
technology (NMR) for the manipulation of the spins, and some rudimentary computations have been performedEI. 
Readout is the main difficulty with this approach, since the field created by a single moment is tiny, and pure states 
Q ' cannot be achieved in the macroscopic samples used. 

, Electron spins also interact weakly with the environment in some circumstance^: relaxation times in excess of 10^ s 
^ ■ have been measured for donor bound states of phosphorus-doped silicon (Si:P)u. The corresponding manipulation 
• I— I . technology (ESR) has also reached a high level of sophistication, but the magnetic moment of the electron exceeds 
' that of the nucleus by 3 orders of magnitude, which presents problems-jof isolation. Readout should be easier than 
^ [ for nuclei, since detection of single electron charges is certainly possibleB, and spin-dependent detection is not far out 
■ - - ' of the reach of current technology. 

Solid-state implementations of QC are particularly attractive because of the possibility of using existing computer 
technology to scale small numbers of qubits up to the 10^ or so that would be needed for nontrivial computations. 
The first paper to propose using the electron spin in a quantum dot subjected to a strong DC magnetic field was that 
of Loss and DiVincenzoO. KaneQ proposed employing the nuclear spin in the Si:P_system as the qubit. A specific 
structure consisting of silicon-germanium (SiGe) layers was proposed by Vrijen et aln. This structure incorporates the 
idea that the factor of an electron can be changed by moving it in a Ge concentration gradient, allowing individual 
electron to be addressed by the external AC field. A different SiGe structure has been proposed by Friesen et alU- 
This structure is designed so that the electron number on the dots, and the coupling between the dots, can be carefully 
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controlled. 

Solid-state implementations must also face the isolation-control dilemma. Decoherence times must exceed the 1 ms 
number in the actual physical structures that are needed for the operation of quantum algorithms. In this paper, 
we examine whether this can be the case for some of the existing proposals based on electron-spin qubits. In the 
process, we hope to learn something about modifications to these structures that can increase r. We shall focus on 
low temperature operation, since, as we shall see, this will probably be necessary in order to obtain sufficiently large 

T. 

We can build on a large body of work, both theoretical and experimental, from the 1950's and 1960's on ESR in 
doped semiconductors. In a series of papers, Feher and coUaboratorsQ investigated the relaxation time for the spin 
of electrons bound on donor sites in lightly doped Si. At sufficiently low temperatures, the relaxation time Ti is 
dominated by single-phonon emission and absorption. In the presence oinSpin-orbit rcoupling (SOC), this can relax 
the spin, causing decoherence. The theory was worked out by Hasegawall3 and Rothllil. 

It must of course be recognized that this spin-lattice relaxation time is not necessarily to be identified with the 
decoherence time. The decoherence time is the shortest time for any process to permanently erase the phase 
information in the wavefunction. This may mean the phase for a single spin, but it also means that the relative 
phases of the wavefunctions of different spins must also be preserved, so that processes that cause mutual decoherence 
must also be taken into account. The actual decoherence time is the minimum of all of these times. A spin relaxation 
time in excess of 1 ms is a necessary, not a sufficient, condition for the viability of a solid-state electron-spin QC 
proposal. 

A QC must have precise input as well as an accurate algorithm. Preparation of the spin state is often proposed 
to be done by thermalization of the spin system at a low temperature. The time to do this actually sets an upper 
limit on the relaxation time of whatever processes thermalize the spins to the lattice. A limit of perhaps 1 — 10 s is 
a reasonable requirement. 

This paper focuses on Ti , the time for relaxation of the longitudinal component of the magnetization, by spin-phonon 
interactions. These processes cause real spin-flip transitions. They occur at random times and thus indubitably 
cause decoherence. In addition to these processes characterized fully by Ti, there are processes which introduce 
random phase changes in the spin wavefunctions. To characterize all such processes bjj-a single "dephasing time" 
T2 will usually not be sufficient for understanding the operation of a multi-qubit systeniE3. Difficulties of definition 
arise, and care must be taken to specify which phase is involved and to what extent it is randomized. For a single 
spin system there is no ambiguity. The 2x2 density matrix pij for the qubit with cylindrical symmetry involves 
only 2 independent parameters pn — P22, and pi2. The time dependence of pn — P22 after a system preparation is 
exponential with decay constant Ti and represents the return of the longitudinal component of the magnetization 
to its equilibrium value. The decay is due to inelastic transitions of the type calculated in this paper. pi2 on the 
other hand, is nonzero only if the preparation of the spin state has a transverse component: Sx{t = 0) ^ 0. The 
decay of this quantity represents the irreversible conversion of this state to an incoherent mixure of "up" and "down" 
states. Again, this is genuine decoherence of the spin state, since the phase information cannaL.be recovered. The 
time dependence of pi2 when the spin is in a strong field has been calculated by Mozyrsky et a^Ej using a Markovian 
approximate master equation. They find that the time scale of the decay due to spin-phonon coupling is very short, 
of the order of the time for a phonon to cross the electron's wavefunction, which is about 10^^"s. But the decay is 
incomplete, with pi2 retaining all but 10~^ of its original value. Their calculation was for Si:P, but a very similar 
result should hold for the dot case. Decoherence at this level is certainly acceptable for quantum computing. These 
authors also pointed out that the decay of the remainder of pi2 is due to the spin-flip processes computed in this 
paper. If T2 is defined as the dominant decay time of the off-diagonal density matrix element, then Ti = T2 for 
spin-phonon processes. 

A quite different source of decoherence is the hyperfine coupling to nuclear spins. The nuclear spins produce an. 
effective random magnetic field on the electrons. Very recent calculations using semiclassical averaging techniquescil 
have obtained a very short relaxation time Ta ~ 1 ns for GaAs-based dot systems in a strong field. This represents the 
decay of the transverse magnetization of an ensemble of dots. This is a dephasing time, but not a decoherence time. 
The electrons spins precess in what is effectively the frozen field of the nuclei. This field is spatially random, and 
the differential precession of the electron spins leads to the magnetization decay. However, this is not an irreversible 
loss of the phase information of the collective wavefunction. Spin echo exeriments are a very beautiful demonstration 
of precisely this point. This " inhomogeneous broadening" presents challenges for the calibration and operation of 
quantum computers, but does not destroy coherence. 

Finally, in any implementation based on electron spin qubits, there will certainly exist small interactions between 
the spins themselves. The dipole-dipole interaction, for one, cannot be avoided, and there may be indirect spin-spin 
interactions mediated by the gates. A recenLpaper suggests that these interactions set the fundamental time scale Tm 
for Si quantum dot implementations of QCcJ. These interactions do produce experimental broadening of ESR lines in 
experiments on bulk systems, and this might be taken as decoherence. In our view, however, these interactions do not 
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destroy the coherence of a state. The system is the set of ah the qubits. During the course of a quantum algorithm 
they are collectively in a pure state (in principle). Any decoherence that destroys the purity of the state comes from 
averaging over the unknown states of the environment. The broadening that comes from dipole-dipole interactions 
comes, in NMR and ESR calculations, from averaging over the states of the system itself, which is not an appropriate 
method for calculating decoherence. The effect of qubit-qubit interactions that cannot be turned off is to complicate 
the quantum algorithm. A quantum algorithm is a unitary transformation that must always include the effect of the 
system Hamiltonian (including dipole-dipole interactions) in addition to external operations. In every case except for 
very simple ones, this algorithm must first be computed, presumably with the help of a classical computer. This step 
in QC may be termed "quantum compilation" . The issue that qubit-qubit interactions raise is not one of decoherence, 
but rather whether the determination of the algorithm, the complation step, becomes prohibitively difficult. This 
could happen for two reasons. One is that the interactions are so poorly known that they cannot be corrected for. It 
seems likely that quantum error correction can resolve this difficulty. A second and more interesting possibility that 
the interactions convert the computation of the algorithm itself into a problem that grows exponentially with the size 
of the system. We regard this as an open question and a deep one, that combines many-body theory with algorithm 
design and error correction. We note that in NME^ implementations the interaction between the qubits also cannot 
be turned off, but it can be canceled by refocusingl3. 

In this paper, our aim is to evaluate the importance of spin-phonon coupling as a source of decoherence in quantum 
dot qubits. Fundamentally, the issue is whether the long relaxation times Ti observed at low temperatures in bulk 
Si:P carry over to SiGc dots proposed for QC. In Sec. II we introduce the structures that we are interested in. In Sec. 
Ill the physics of the spin-phonon relaxation mechanism is described qualitatively. In Sec. IV calculations and results 
for pure Si systems are given as a function of the critical design and operating parameters. In Sec. ^ are found the 
correponding results for structures containing Ge. Sec. VI is the conclusion. 



II. THE STRAINED SILICON QUANTUM WELL 



Si-Ge heterostructures are utilized widely in the digital electronics industry, and presently have the shortest switch- 
ing times of any device. One reason for their success lies in the ability to engineer structures of near perfect purity, 
with control over thicknesses and interfaces that approaches atomic precision - a technological tour-de-force. An 
equally key achievement has been the harnessing of strain as a tool to control band offsets in heterostructure device*. 
This paper presents calculations of spin relaxation for real SiGe structures such as those proposed by Vrijen et al.a 
and Friesen et aln. Accordingly, we have calculated the electron wavefunctions ui quantum wells, which is needed as 
input for these calculations. Details of these calculations are presented in Rcf.Ll, and will not be repeated here. In 
this section we only describe those aspects of the calculations that are germane to spin relaxation. 

Quantum wells are constructed by sandwiching a very thin layer of one material between two others. Electrons 
can be confined in the quantum well layer when the conduction band offsets produce a potential well. The key to 
this technology is therefore to understand the band structure of the various layers. In this section we will consider 
a particular class of wells formed of pure Si, sandwiched between barrier layers of SiGe. We will find that this is 
optimal from the standpoint of spin coherence. Metallic gates or impurities create zero-dimensional bound states that 
define a quantum dot. We first review briefly effective mass theory for dots in pure, unstrained Si, then unstrained 
SiGe, and finally strained Si. 

In pure Si, the A conduction band minima occur near the symmetry points X , in the directions {001}. In a perfect 
Si crystal these minima are six-fold degenerate, the valleys being equivalent. In the dot the electron feels a potential 
Vg{r) in addition to the atomic potential, which lifts the degeneracy, though the splittings are not large. The spatial 
variation of Vg{r) is on length scales generally much longer than the lattice spacing. For the moment, we shall assume 
that the electron is in the ground state of Vg{r) and ignore mixing with any excited states. To the extent that the 
scale of variation of Vg (r) is much longer than the atomic spacing there are six pearly degenerate ground states. This 
is referred to as the "valley degeneracy". The wavefunctions can be written asllZl 
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a>„(f) = ^a(^)F,(r>,(r). (1) 

Here, (f>j is a Bloch function of the form 

cl^,{f)=u,(r)^''^-~, (2) 

where kj are the six A minima {+kQX, — fegi, +koy, —k^y, -\-kQZ, —k^z} (we shall always use this ordering), and Uj{r) 
are periodic functions with the same periodicity as the crystal potential Vp{r). The F±z are envelope functions that 
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satisfy the Schrodinger-like equation 

-2;^ a;^ - 2;^ + a?) + ^ - <y±.«, (3) 

and are independentfy normalized to unity, similar to wavefunctions. Analogous equations can be given for the ±x 
and ±y minima. We see that = F^^, Fy = F-y, and F:^ = F-z, so only three independent envelope functions 
must be computed, mi and rrit are the longitudinal and transverse effective masses associated with the anisotropic 

conduction band valleys. E^^^ is the A conduction band edge at h^. The splitting of the degeneracy comes from 

(7) 

corrections to this envelope fuction a ppr oximation. Different choices of the constants aii determine the six states. 



Their values will be discussed in Sec. III. This formalism is a good approximation for both dot and impurity bound 
states, as the valley splitings are much smaller than the energy scales in 

Germanium is completely miscible in Si, forming a random alloy. For variable Ge content x, Sii_a;Gea; exhibits 
materials properties that vary gradually over the composition range. The alloy lattice constant, ao{x), follows a linear 
interpolation between pure Si and Ge, known as Vegard's law, quite accurately for all x: ao{x) = (1 — x)asi +a;aGcc3. 
Electronic properties show an abrupt change in behavior near x ~ 0.85, where the Si-like A minima cross over to 
four-fold degenerate, Ge-like, L minima. In this work we focus on the range x < 0.5, which is strictly Si-like, though 
we will have some remarks below on Ge-rich structures. Throughout this range, properties such as effective mass and 
the dielectric constant vary only slightly from pure Si values. For our calculations, the most important parameter is 
the conduction band edge, E^^\x), which remains 6-fold degenerate in the range x < 0.5. The theory of the variation 
of E'-^Kith X is not-germane to the present work, and we simply quote the empirical result, linear in x, which is 
consistent with Ref.EJ: 

AE^'^^x) = E^'^^x) - i;(^)(0) ~ 0.23a;(eV). (4) 

(We note, however, that Ref.lli suggests a slope for E^^^x) of opposite sign.) The relatively weak variations of the 
effective mass and the dielectric constant will be ignored here. 

We consider thin Si wells in which the Si layer grows pseudomorphically. The in-plane lattice constant, a||, must 
be the same for all the layers, causing a tetragonal distortion in the strained layer(s). Here we consider the case of 
strained Si grown on the (001) surface of relaxed Sii_2,Gej;. The in-plane Si lattice constant depends on x as 

a|| (x) = (1 - a;)asi + xacc- (5) 

Since aco > isii the Si is under tensile strain in the plane. Hence, the out-of-plane Si lattice constant, aj_, is reduced 
according to continuum elastic theory: 



aj_{x) = asi 



^ _ ^ ci2 a\\ix) - asi 
Cll asi 



(6) 



where cn and C12 are elastic constants for pure Si. 

Strain produces shifts of the A band proportional to the strain variables 

a||(.x)-asi a±{x)-asi 

e|l(a;) = — and e_L(a:) — . (7) 

asi asi 

with proportionality constants called the dilational and uniaxial deformation potentials, S^'^'' and respectively. 
Because of the anisotropic nature of the strain, the two z minima are shifted down relative to the x and y minima, 
restating in a splitting of the A conduction band. The net shifts with respect to the unstrained Si A band are given 
byM 

AE^^-\x) = (4'^) + is('^))(2e||(x)+£^(x)) + ^S('^)(£^(x)-£||(x)), (8) 
Ai?(^ii)(a;) = (4'^) + is('^))(2e||(x)+£^(x))-isf)(£^(x)--£||(x)). (9) 

The first terms in Eqs. (H) and (^) are hydrostatic strain terms, which lower the conduction edge compared to 
unstrained Si. The second terms in Eqs. (|^) and (^ produce the splitting, associated with uniaxial strain. To perform 
our calculations, we use the materials parameters given in Table |. However we note that the deformation potentials, 
particularly S^'^'' , are very difficult to measure experimentally. Considerable disagreement exists in the literature as 
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to the value and even the sign of 'c3. The value given in Table,| was-j-eported (but not endorsed!) in RefJi, and 
provides energy band variations in general agreement with Refs.tSi andta. We arrive at the following strain-induced 
shifts of the conduction band edge for pure Si: 



^^(Ai) ^ _o.86a; and AE^'^^^^ = -0.16a 



(10) 



The corresponding shift in the relaxed barrier layers, due to the presence of Ge, was given in Eq. (1). Together, these 
results describe the conduction band offsets for the quantum well that are used in our simulations. 

We now apply our results to two specific quantum well designs of intw'est for quantum computing. Design (1), 
shown in the inset of Fig. |l], is a version of that proposed by Vrijen et al.u, in which electrons arc trapped on donor 
ions (usually P), implanted in a semiconductor matrix. In that work, the quantum well is split into Ge-rich and 
Si-rich regions to facilitate single qubit operations. For simplicity, we consider here a uniform quantum well, formed 
of pure Si, with a single dopant ion located ai—the center of the well. In such a device, jSingle qubit operations can 
be accomplished using a coded qubit schemeEj. Design (2), proposed by Friesen et alu, is shown in Fig. The 
confinement potential for the electrons is much softer than in design (1). Electrons are trapped vertically by the 
quantum well, and laterally by the electrostatic potential arising from lithographically patterned, metallic top-gates. 
Additionally, the quantum dot is tunnel- coupled to a degenerate doped back-gate. The dimensions for both designs 
are given in the figures. 

The wave function of the bound electron is computed in the envelope function formalism, Eq. (^. Couplings 
between the different valleys are introduced through the perturbation theory described in Sec. |V|. This procedure 
provides the specific values of for the ground state, which we use in our calculations. 

The abrupt conduction band offsets are handled by matching the ground state wave function ^g{r) and dz^g{r) 
at the interfaces. (Remember that we have equated effective masses on both sides of interfaces.) Due to linear 
independence of the Bloch functions, the boundary conditions do not cause mixing of the envelope functions. Solutions 
of Eq. (H) and the analogous equation are obtained, using commercial 3D finite element software. 

As will be seen in Sec. ^ the key quantity for the computation of Ti is which describes the probability for the 
bound electron to be on a Ge atom. Ge is associated with reduced coherence times, by virtue of its large spin-orbit 
coupling. Referring to Eq. (1), one deduces that this probability may be expressed as 



4(4-))2/' dVFj(r) + 2(4^))2 /■ d^F^ir) 



(11) 



where fit, is the volume outside the quantum well, if the well is pure Si. The subscript g refers to the ground state. 
The term in square brackets reflects the probability of finding the bound electron in a barrier region, while x gives 
the probability that the electron is on a Ge site. 

Fig. |l| shows the results of our calculations for in design (1), as a function of the Ge concentration, x, in the 
Sii_a;Gea; barriers. For x > 0.02, decreases with x for two reasons. First, as x increases, the conduction band offset 
at the quantum well also increases, allowing less of the wavefunction to penetrate the barrier. Second, the spatial 
extent the electron in the z direction is greater for F^ than Fz, because of the anisotropic effective mass . However, 
less of Fx is mixed into the wavefunction for large x, since ax becomes very small. Eor x < 0.02, drops quickly to 
zero, due to the absence of Ge in the barriers. In the actual design of Vrijen et al.p there is Ge in the active layer. 
To give an idea of the effect of this, we include an equivalent value of for such a structure. 

Fig. ^ shows results of the calculation for design (2), as a function of the quantum well thickness, z. To perform 
the calculations, we have considered a fixed Ge concentration, x — 0.05, and taken the limit of large strain, so that 



{±x) 



and ckg*^'' ~ l/\/2 for the ground state. As z increases, less of the wavefunction penetrates the 



barrier regions, causing / to decrease. 



III. SPIN RELAXATION DUE TO COUPLING TO PHONONS 



In this section we give the method for calculating Ti , the spin- flip time of a spia qubit in tiije ground orbital state 
due to emission or absorption of a phonon, following the logic used by HasegawaEiJ and Rothtll for bulk Si. 

Consider a single impurity with a unit positive charge, such as a phosphorus atom, at the origin. In the absence 
of central cell corrections, there is a 12-fold degenerate ground state, including spin. This valley degeneracy of the 
ground state is reduced to two by these corrections, and the splitting between the 2-fold spin-degenerate ground state 
and the higher states is of order Ai? ~ 10 raeV. We shall discuss the detailed linear combinations (a- values) of the 
states below, as the coefficients giving the various valley amplitudes play an important role in the calculation of matrix 
elements. These twelve states may all be thought of as hydrogenic Is states. The splitting of Is and 2s is about 
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30 meV , larger than the 10 meV valley splittings. Let us now split the 2-fold degenerate ground state by applying a 
DC magnetic field in the z-direction. The transition rates between these states are denoted by and The 
relaxation time Ti is defined by 1/Ti = W^^l + M^x|. 

The transitions are caused by phonons, but there are important approximate symmetries that suppress these 
transitions. These are: (1) spin rotation symmetry, meaning that the electron spin cannot be fiipped by a phonon; 
this symmetry is broken by SOC; (2) time-reversal symmetry, meaning that one state cannot be changed into its 
time-reversed partner by emission or absorption of a phonon; this symmetry is broken by the external magnetic field; 
(3) point group symmetries; these are partially broken when strain is applied. 

The spin rotation symmetry would rule out phonon-mediated transitions between the two states entirely if there 
were no spin-orbit couphng (SOC). This means that the effects of SOC on the wavefunctions, even though these 
effects are small in relatively \ow-Z Si, must be taken into account. When we refer to a state as f or J,, these symbols 
must be taken to refer to the majority spin content of the state, not to a pure spin state. Transition rates are roughly 
proportional to {g — 2)^, [more precisely {gi — gt)"^, where gt{gi) is the transverse (longitudinal) g-factor, see below for 
definitions] . 

The time-reversal symmetry implies that transitions cannot take place directly between Kramers-degenerate states 
even in the presence of spin-orbit coupling. The direct phonon-mediated transitions between the two states of interest 
to us are strongly suppressed by this approximate symmetry. It is broken only by the external field H. The fastest 
processes then involve virtual excitation to higher-energy states that are mixed into the ground state by H. Hence 
l/Ti involves a factor {^ibH/AE) . There is an additional factor of from the phonon density of states, giving an 
overall rate l/Ti ^ in the limit of small H. 

The point group symmetry is reduced from cubic to tetragonal under strain. This has complex effects that we will 
explain below. 

Before giving actual calculations, we summarize those differences between the electrons in donor impurity states 
and in an artificial dot that affect Ti. The most obvious is the single-particle potential that binds the electron. The 
gate potential is much smoother than the hydrogenic potential of the impurity. This implies that the corrections 
to the effective mass approximation are much weaker, and AE will be much reduced. It is difficult to compute the 
energy splittings precisely, but considerations based on the method of Sham and NakayamacJ give splittings in the 
range ^ 0.05 — 0.1 meV in the structure of Friesen et al.El. This increases the relaxation rate. (In fact a naive estimate 
of the enhancement is a factor of 400.) On the other hand, the structures we consider have strong lattice strain. This 
partly lifts the valley degeneracy and also reduces the matrix elements, which decreases the rate. Another aspect of 
some of the proposed designs is the presence of Ge with its much stronger SOC. This will act to decrease the spin 
relaxation time. 



IV. PURE SI QUANTUM DOTS 

We consider first the case of pure Si under uniaxial strain. The ingredients of the calculation are as follows. 

From Sec. II we have the solutions to the Schrodinger equation [Ho + Vg(r)]$„o(r) = i?„$„o(^)- '^o is the 
unperturbed crystal Hamiltonian without SOC and it has the full space group symmetry. Vg{r) is the gate and/or 
impurity potential. 

To calculate Ti, we must also include SOC, which we treat as a perturbation: Hsoc = ^SiJ^r^r ' ^r- The 
resulting states $n(0 two- fold degenerate because of time- reversal symmetry. Let us denote these states as 
^nii'P) and $„|(r). They are not eigenstates of spin, so the arrows denote pseudo-spin. We may define the pseudo- 
up state as that which evolves from the spin-up state as spin-orbit coupling is turned on adiabatically, and similarly 
for the pseudo-down state. Because of valley and pseudospin degeneracy, there are 2 ground states ^gs{r) and 10 
excited states ^rs{r)- The 12-fold degeneracy in the effective mass approximation is broken by central-cell corrections 
in the impurity case and smaller corrections in the quantum dot. 

There is also the Zeeman Hamiltonian of the external field Hz = I-i-bB ■ (L + 25). In the field <I>„|(r) and $„|(r) 
are no longer degenerate. Note that the energy splitting may depend on the direction of B. 

Finally, we have the electron- phonon coupling Hamiltonian Hep- A phonon represents a time-dependent perturba- 
tion. This will create transitions whose rate is given by the Fermi golden rule. We are interested in the transitions 
between $c,|(r) and $g|(r). However, {^gi{r)\Hep\^gi{r)) = in the absence of the external field. Thus we need 
to calculate in next order in perturbation theory using an effective Hamiltonian 

^' = p \ {[Hz\^rs{r^) {^rs{r)\Hep] + [Hep\<^rs{7^) (f) |Hz ] } (12) 



Here r runs over the excited states, r = 2, .., 6. s =t) i- 
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The relaxation time is given by 1/Ti = W^^i + Wx|, where Wi^ = {27r/h) x J^qX I (*st(^I^' -^(^gi - 

Eg-\ — hw^x) [1 + n{uj^x)\ is the rate for transitions from the higher-energy (pseudo-spin down) state to the lower- 
energy (pseudo-spin up) state and W^i = (27r/?i) Yl,qx I ( ^si(^) l^'l *''st(^))I^ '^(-^si ~ ^sT ~ ^qx) n{uj^x) is the rate 
for transitions from the lower-energy (pseudo-spin up) state to the higher-energy (pseudo-spin down) state where the 
sum is over phonon modes qX with energies uj^x ■ A thermodynamic average over the lattice states has been taken. 
It yields the Bose occupation factors n{ujq'\) for the phonons. 

The matrix elements of Ti' are computed as follows. 

The expectation value of the external field Hz can be written as 
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.1=1 



= ^ibB ■ D„„/ • as,s' (13) 

and the tensor g(') is the effective g factor at the ith valley. This equation defines the tensor D„„' that characterizes 
the coupling of the various states by the external field. The principal axes of the g tensor are the same as that of the 
effective mass tensor at the ith valley. It has the form 

gi 0\ / 5t \ / 5t 

g(±-) = I , =05/0 ,g^^^^ = | . (14) 

gt J \0 gt J V 5/ 

There are only two independent constants. 
A simple example of the diagonal part of the D tensor is that for the ground state of an impurity in the unstrained 

lattice when the central cell corrections are included. Then we have a^^ = I/VG and Dgg is proportional to the unit 
matrix: 

{^gs{r^\H,\^gs,ir)) = ggtXBB ■ dss' , (15) 

with 

2 1 

9g = 3.9* + ^91- (16) 
The matrix elements between ground and excited states have the form 

($g,(f)|iJ2|$,,,(f)) ^ g' ^1BB ■ Bgr ■ <7,s' , (17) 

with 

9' = l{9i~9t) (18) 

and the tensor Dr is defined by 

6 

Dg, = 3^a«a(')W« (19) 

i=l 

where fc*^*^ is the "local" anisotropy axis. If the original g were isotropic, then g' = and there would be no coupling 
between different states and no spin relaxation. 

If the lattice is strained, then the a coefficients become strain-dependent and the general expression for D from Eq. 
[Tsl must be used. Uniaxial strain lifts the degeneracy of the valleys. We include this effect in the Hamiltonian and it 
determines the proper combinations of the al*' defined in Eq. ^ These then feed into D^j.- As a function of strain 
the ckn^'* for the ground state cross over from the completely symmetric combination a^'' = to the combination 

Qfg^^'' — = 0, = I/V2 in the limit of large strain.a 
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The phonons involved are just the acoustic ones, one longitudinal and two transverse - these are the only ones with 
low enough energy to play a role in relaxing the spins. The matrix elements of the electron-phonon interaction are 
only nonzero within one valley and for a single phonon mode they are conventionally parametrized as: 



« (f)|i7(f = ^b,^,eM ■ (Sdl + S„fc«fc«) • q+ h.c. 



(20) 



near the ith valley, where q = k — k' and e\ is the polarization vector. destroys a phonon with wavevector q and 
polarization A. Once again, we see that the interaction can be characterized by just 2 parameters, in this case and 
S„, as already defined in Eq. ^ Performing the integration over the envelope function at wavector q now gives 



$,,(f)|ij(f)|$. 



[r]) ^ = A{q) {iex{q) ■ T>gr -qb^ix- ie*x{q) ■ Dg^ • q b*^x) 



(21) 
(22) 



where 



A«(q) = Y^F^^>{k + q)F^'\k) = J d^rF'^{r)e'' 



(23) 



Thus the electron-phonon interaction involves a form factor for the bound states. Since F is normalized, we have 
A^*-' ((j) « 1 when the wavelength of the phonon is much longer than the spatial extent a* of the bound state: qa* << 1. 
The calculations of II indicate that this is the case. It is also independent of (i). 

In the golden rule calculation, the energy dcnoninator {Eg — Er) will suppress contributions from the excited 
states of Vg. Thus we will keep only states that are split off from the ground state by corrections to the effective 
mass approximation. This approach works very well in Si:P and should be even better for the quantum dot. 

This produces the golden-rule transition rate 



2n 
li 



A'^iq) 5{Eg^ - Egi - hw^x) {aq\a}x) 

q\ 

2 



E 



B ■ Tigr ■ (T|X ^^(q) ■ Dgr • qSs 



Eg — Er 



(24) 



We approximate the phonon dispersion as lo^\ = vxq. Setting A — \, performing the integral over the magnitude 
of (f, and repeating the calculation for IV||, wc obtain a total spin relaxation time 



A=l >• "'0 

Here {6, (/>) are the polar axes of the direction of B measured from the [100] direction of the crystal, p is the mass 
density. This is our basic result, in a form very similar to that given by HasegawaEiJ. 

In Fig. ^ the effect of strain on the lowest energy levels is shown. Here the zero-strain splittings are those of a P 
impurity. A similar plot was given in Ref.cl, but at that time the correct zero-strain splittings were not known. In 
the dot, the energy splittings for the unstrained case would be two to three orders of magnitude smaller (see below 
for a discussion). 

In Fig. ^, we show the effect of strain on the relaxation time of an electron in a P impurity potential. There are 
two effects: the overall increase of the energy denominators and the change of the ground state to a less symmetric 
valley weighting. This leads to a non-monotonic dependence of Ti on strain, but in the region of interest, the effect 
of strain is to greatly increase Ti, since most proposed structures have s < — 1. At large strain, only one energy 
denominator remains small, that between the ground state, symmetric in the ±z valleys, and the first excited state, 
antisymmetric in the ±z valleys. The overlap matrix D is very small between these two state. This reduction of the 
matrix element is the dominant effect. 



t> 

n 

r=2 



B( 



D 



gr 



exiq)--D 



Eg — Er 



(25) 
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Also of interest is the dependence of Ti on the angle of the external field, since this may serve as a diagnostic tool 
in experiments to verify that the relaxation process is really due to spin-phonon coupling. As a function of strain, 
this dependence becomes highly anisotropic, as seen in Fig. ^. This is due to the elimination of all but the zLz valleys 
from the problem at high strain. We note that the limiting value Ti — > oo when the,external field is along a crystal 
axis is cut off by intcrvalley scattering effects not included in the present calculationQ. 

The change in the confining potential reduces the corrections to the effective mass approximation, as discussed in 
Sec. [1] which in turn reduces the energy denominators, leading to a decrease in Ti.jJji the Friesen et al. structure, 
we have estimated the splitttings Er — Eg using the method of Sham and Nakayamao, and they range from 0.05 to 
0.1 meV, depending on the gate voltages that produce the potential. This may be incorporated into the calculation 
of the a'n' by introducing a variable coupling that mixes the ±x and iy valley wavefunctions with the ±z valley 
wavef unctions. For a precise definition of Ac, see Ref.a. Approximately, however, Ac ~ 0.2 [Er — Eg) < 10~^eV. 
The results of the calculations are shown in Fig. ^. The graphs show the absolutely crucial role that strain plays in 
the determination of Ti. Because of the small energy denominators in the dot, the rate is extremely fast at small Ac 
for the unstrained case: 1/Ti ^ ^J^- 

We also point out that in stuctures where Ti appears to be too long for efficient preparation of the spins (as discussed 
in Sec. I) the deficiency can be made up by increasing the temperature T. As seen in Fig. |^, Ti decreases very rapidly 
as T is increased. Actually, this calculation even underestimates the decrease, since multiphonon processes begin to 
contribute at about 3K. However, it must be borne in mind that the temperature must be small enough that the 
initial state is essentially pure (all spins up, for example). A more practical method of decreasing Ti if desired, would 
be to increase the magnetic field. 



V. STRUCTURES CONTAINING GE 

The effect of alloying with Ge is to increase the SOC and hence to increase — 2|. First principles calculations 
of Ge impurities in Si have shown that there is little effect on the statesj-uear the bottom of the conduction band, 
though this is not necessarily the case for higher energy states in the bandEj. This is in accord with the isoelectronic 
character of the atoms in the alloy. It is then reasonable to employ the virtual crystal approximation (VGA). The 
approximation should be quantitatively accurate for small, (say < 10%) concentrations of Ge, but may be taken as a 
good qualitative guide also to higher concentrations. The Bloch functions in the absence of SOG satisfy 

^o*„£,=i?A£„ (26) 

where 

^ 2 3 

= ^^exp{ik ■ rji) ^ ai,n{nk) 4i„is{r - r^i) (27) 



* r = ^ 

^ j = l 1=1 m=0 

where Nc is the number of unit cells, j labels the unit cells, / labels the two positions in the unit cell, m labels the 4 
atomic states, (j)ms are the atomic orbitals, and the coefficients aim{nk) give the proper linear combination of atomic 
orbitals (LGAO's) for the state at momentum k and band n. The normalization condition is 

Y,\aim{nk)\^ ^1. (28) 

Im 

In the VGA, the aim{nk) are the same for the Si and Ge sites. The SOC Hamiltonian for the ahoy is written as 

i?soc = E^^^«-% (29) 

R 

where A^^ is the SOG strength for Si (Ge) when is a Si (Ge) site. The energy shift of an electron in the external 
field B = (Ba;, By, Bz) can be written as 

AE^^^(B) = -^iB J2 9^{nk)B,s, (30) 



i—x.y,z 
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TABLE I: Materials parameters used in this work. 



where 



Parameter Value 


ft 


Pf- [# 1 


asi 5.43 A 


Q 




flGc 5.66 A 


1£ 




cii 1.675 


It 




ci2 0.650 


IS 




Si^^ 9.29 eV 


If 




h'^> -10.7 eV 


21 





gi{nk) = 2 - 2fi.2 



E 



air. 



(nk) aira'{n'k) {E^j: 



(31) 



and Eimm' is the completely antisymmetric symbol with £123 = 1. 

Our approximation now consists in regarding the dependence on the Ge concentration as coming only in the 
term {^j^^j^/Nc) ■ Hence any component of the g-tensor is proportional to a weighted average of Xsi and Xce- 
In particular, consider the conduction band minimum at k = (0,0, fco) in the compound Sii-^^Gca;. Then we have 
gz{x) — 2 = [gz{Q) — 2] [(1 — x) + x{\Ge/Xsi)\ , where 5z(0) is the value for pure silicon. The numbers, known from 
atomic physics, are Xce/^Si — 0.295/0.044 = 6.71, so gz{x) — 2 — [gz{0) — 2] [(1 — x) + 6.71a:] . Similarly, we have 
gy{x) — 2 = [gy{0) — 2] [(1 — x) + 6.71x] , etc. Only the overall scale of g changes, not the anisotropy in the VGA. 

The presence of Ge in the lattice breaks the translation symmetry, an effect that is neglected in the VGA. Taking 
this into account would lead to a (7-factor that is averaged over momentum space rather than one that is evaluated 
only at the Si conduction band minima. If we take the L-minimum of pure Ge as an example, we would expect that 
the corrections to Iff — 2| would be greater if there is some averaging over momentum space. Thus a calculation going 
beyond the VGA would give corrections that would reduce Ti. 

We show the results for differing Ge concentrations in the active region of the well in Fig. ^. In the Ge-poor regime 
X < 0.85, the relaxation rate 1/Ti decreases fairly slowly. Thus designs with a Si-rich well will not suffer from very 
fast spin relaxation on the minority Ge sites. 

By contrast, if we cross the critical point around x = 0.85 where the conduction band minima switch positions, 
then there is a very rapid decrease of Ti. This is due to the much greater g- factor anisotropy of the [111] minima: 
g' w —0.4, (as compared to g' < 0.001 in Si). The results for the Ge-rich region are much simpler than in the Si-rich 
region. Uniaxial strain does not affect the relative valley energies or the matrix elements parametrized by D. The 
only effect of going to the dot from the impurity is to reduce corrections to the effective mass approximation, which 
will strongly decrease the relaxation time. Indeed, if we take the impurity result for Ge:P from Ref.E3, which is 
Ti = 2.3 X 10"'^ s and estimate the decrease in AE as about a factor of 100, then we obtain Ti - 10"^ to lO-'^s, 
which is vastly shorter than that of Si-rich structures. The physics is the same as that which governs the divergence 
of the rate in Fig. ^ for the unstrained case, namely the much smaller energy denominators. 

To make contact with Sec. II note that if there is Ge in the barrier regions, as is the case in most SiGe designs, then 
we will have an envelope function that weights Si and Ge differently, then x in the VGA formulas must be replaced 
by — X^fl l/flPj where is the amplitude that the electron is on a Ge site and the sum runs only over Ge 

sites, as described in Sec. II. This is the product of the Ge concentration in the layer, times the probability that the 
electron is on the layer, summed over layers. Referring to Figs. |^ and ^, we see that is always less than 10^"^ for 
structures that have pure Si active layers and Ge in the barrier layers. The effect on Ti is small so these structures 
are "phonon-safe" . 



VI. CONCLUSION 



Decoherence due to spin relaxation by emission and absorption of single phonons does not pose a major obstacle to 
quantum computation in quantum dot SiGe heterostructures using spin qubits, if these qubits are properly designed. 
This spin relaxation mechanism is the dominant one in the case of electrons in donor bound states, such as in Si:P. 

In order that this source of decoherence be kept under control, however, certain conditions must be satisfied in the 
design of the structure. Unstrained SiGe alloys generally have relatively short spin relaxation times due to valley 
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degeneracy. This degeneracy produces strong spin mixing in the Kramers-degenerate ground states of the quantum 
dot once the field is turned on. This in turn increases the electron-phonon matrix element that connects the two states 
and decreases the relaxation time. This effect, fundamentally due to spin-orbit coupling, is stronger in the quantum 
dot than in the impurity state. The corrections to effective-mass theory, stronger for the impurity potential, tend to 
lift the degeneracy. 

The crucial role that strain plays, in Si-rich structures, is to lift the valley degeneracy. In this respect compressive 
uniaxial strain is the best, since the residual degeneracy is reduced to two valleys. The wavefunction for the ground 
state is then symmetric in the valley index, and the first excited state, the only one with a small energy denominator, 
is antisymmetric in this index. This fact strongly suppresses the spin mixing effect of the external field. Thus a 
workable design should include uniaxial strain that exceeds a critical value. This value is small enough that there will 
be sufficient strain in most practical structures in which the active layer is Si rich and sandwiched between relatively 
Ge-rich layers. The active layer must be thin enough to avoid dislocation-mediated relaxation. 

Spin-orbit coupling is stronger in Ge than in Si due to the higher atomic number. This suggests that Si-rich 
structures are to be preferred, and the calculations bear this out. The relaxation time decreases with increasing x 
in a Sii-ajGca; alloy. In fact, if x increases beyond x — 0.85, decoherence becomes very strong, and quantum dots in 
these Ge-rich structures, such as those proposed by Vrijen et alH, may well run into difficulties for this reason. A 
key point for these structures is that uniaxial stress does not lift the valley degeneracy, since the valleys have moved 
off the crystal axes. Strain-induced suppression of spin relaxation cannot occur. 

Beyond design issues, fabrication quality is also important. Even a small concentration of magnetic impurities will 
negatively impact the spin relaxation. In modern semiconductor technology, however, concentrations of magnetic 
impurities much less than 0.1 ppm are routine. The impact of lattice imperfections is less clear. These will generally 
act to lower the symmetry of the system, and to lessen the accuracy of the effective-mass approximation. As we have 
seen in the context of the impurity calculations, these effects generally increase the spin relaxation time because of 
reduced degeneracy. On the other hand, localized bound states can form at such imperfections. If this results in 
free electron spins, it could have seriously negative effects on the relaxation time. Similar effects would result from 
lattice defects that produce localized phonon modes. 

Finally, it appears that quantum dot designs will need rather low temperatures in order to ooerate. The impurity 
experiments clearly indicate that multi-phonon relaxation increase rapidly with temperaturecj, with a crossover to 
this regime at about 3 K. We speculate that this crossover temperature is pushed up when the effects of strain are 
included, but we have not yet done the requisite calculations. 

We would like to acknowledge useful discussions with our collaborators at the University of Wisconsin-Madison 
and with Xuedong Hu. This work was supported by the National Science Foundation under the Materials Theory 
Program, Grant No. DMR-0G81039, the Materials Research Science and Engineering Center Program, Grant No. 
DMR-0079983, and the Army Research Office. 
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FIG. 1: Probability, f^, for finding a donor-bound electron at a Ge site, as a function of Ge content, x. The simulated structure, 
design (1), is shown in the inset. A strained Si quantum well of thickness 6 nm is sandwiched between relaxed Sii-i^Gea; barrier 
regions of thickness 20 nm. The electron is bound to a ion embedded in the center of the quantum well. 




Well thickness, z (nm) 

FIG. 2: Probability, for finding an electrostatically bound electron at a Ge site, as a function of quantum well thickness, 

z. The inset shows the heterostructurc layers for design (2), beginning at bottom: a thick, doped semiconductor back-gate, 
a relaxed Sii_a;Gea; barrier layer, a strained Si quantum well, a thick, relaxed Sii-xGea: barrier layer, and lithographically 
patterned metallic top>-gates. The distance between back- and top-gates is held fixed at 40 nm, while the quantum well, of 
variable thickness, z, is centered 15 nm above the back-gate. 
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FIG. 3: Energy of the six Is-like donor levels of an electron bound by a P impurity in silicon (with respect to the energy 
center of gravity) versus the strain parameter s with an uniaxial stress applied in the [100] direction. (It is important to note 
for clarity that what we call compression relative to the [100] or growth direction is equivalent to tension in the plane of the 
quantum well. This label convention differs across fields.) For reference, s = —3 corresponds to the compressive (s negative) 
strain caused in a pure silicon layer by a Sio.8Geo.2 sublayer. The energies are expressed in eV and the numbers in parenthesis 
indicate the degeneracy of the level. The Ai (ground state, solid line) and Ti level wave functions are mixed by strain (causing 
valley population intermixing) and the relaxation rate is proportional to 1/Ai?. Thus,[-ivith all else equal, an increase in strain 
causes the riflsixation rate to decrease. For a detailed analysis see Wilson and Feheiu, where s here is their s' times 100 or 
Koiller et aZ.Ea whose notation is s = 100(6xAc/Hc). 
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FIG. 4: Relaxation rates of a P impurity bound electron in [100] uniaxially strained silicon versus strain parameter s for a 
temperature of 3 Kelvin and a magnetic field H — _ff (cos S cos 0, cos sin of 1 Tesla in the [111] {9 = (f> = 7r/4) and [110] 
{9 — it/4:, = 0) directions respectively, s = —3 corresponds to the strain caused in the pure silicon layer by a Sio.8Geo.2 
sublayer. 
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FIG. 5: Dependence of the relaxation rate on magnetic field direction with H = H {cos 9 cos (j), cos 0, sin 9) for (a) unstrained Si 
and (b) compressively strained Si (corresponding to a Sio.8Geo.2 sublayer). The magnetic field is set at 1 T and the temperature 
is 3 K. 
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FIG. 6: Relaxation rate of a bound electron in silicon. The solid line is for fixed [100] uniaxial strain s — —3, corresponding to 
the strain caused in the pure silicon layer by a Sio.8Geo.2 sublayer. The dashed line is the unstrained case. The rate is plotted 
for varying the intervalley coupling constant Ac. Ac corresponds to l/6th the energy splitting between the singlet (ground 
state) and triplet (next excited state) valley energies in unstrained Si. As a result, it controls the mixing between the ±j: and 
±1/ valley wavefunctions with the ±z valley wavefunctions. The magnetic field is set at 1 T along the [111] direction and the 
temperature is 3 K. 
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FIG. 7: Relaxation rates of a P impurity bound electron in [100] uniaxially strained silicon (s = —3, corresponding to the 
strain caused in the pure silicon layer by a Sio.sGeo.2 sublayer) versus (a) temperature with the magnetic field set at 1 T and 
(b) magnetic field strength with the temperature set at 3 K. 
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FIG. 8: Relaxation rate of a P impurity bound electron in [100] uniaxially strained silicon versus concentration of germanium 
within unstrained (left vertical axis) and compressively strained (right vertical axis, corresponding to a Sio.8Geo.2 sublayer) 
pure silicon. The magnetic field is set at 1 T along the [111] direction and the temperature is 3 K. Here we assume that the 
addition of germanium does not affect the strain within the bulk and only acts to increase the g-factor and spin orbit coupling. 
We expect this approximation to hold for small concentrations of germanium and to break down with increasing concentration 
as the conduction band valley minima switch from the X type of silicon to the L type of germanium. 



